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We outline the qualitatively different physics behind charging-induced current asymmetries in 
molecular conductors operating in the weakly interacting self-consistent field (SCF) and the strongly 
interacting Coulomb Blockade (CB) regimes. A conductance asymmetry arises in SCF because of 
the unequal mean-field potentials that shift a closed-shell conducting level differently for positive 
and negative bias. A very different current asymmetry arises for CB due to the unequal number 
of open-shell excitation channels at opposite bias voltages. The CB regime, dominated by single 
charge effects, typically requires a computationally demanding many-electron or Fock space descrip- 
tion. However, our analysis of molecular Coulomb Blockade measurements reveals that many novel 
signatures can be explained using a simpler orthodox model that involves an incoherent sum of Fock 
space excitations and hence treats the molecule as a metallic dot or an island. This also reduces 
the complexity of the Fock space description by just including various charge configurations only, 
thus partially underscoring the importance of electronic structure, while retaining the essence of 
the single charge nature of the transport process. We finally point out, however, that the inclusion 
of electronic structure and hence well-resolved Fock space excitations is crucial in some notable 
examples. 
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I. INTRODUCTION 

Ever since its inception [l| , molecular rectification con- 
tinues to be of great practical interest. While rectification 
could arise from asymmetries in the intrinsic molecular 
structure or vacuum barriers at the ends, there are mul- 
tiple experiments 0, [3, Q that exhibit pronounced asym- 
metries in current voltage (I-V) or conductance voltage 
(G-V) characteristics due to unequal coupling with con- 
tacts. Fig. [T] shows that the nature of contact-induced 
asymmetry is qualitatively different depending on the 
nature of the molecule-contact bonding. For molecules 
strongly coupled with the contacts with adiabatic charge 
addition, equal current plateaus are reached over unequal 
SCF voltage widths (Fig^ [T^) , leading to prominent con- 
ductance asymmetries [2| . The origin of this asymmetry 
is the different average charging energies that generate 
unequal mean-field potentials for opposite bias voltages 
[1]. Reducing the contact-molecular coupling drives the 
system into CB, where the intermediate open-shell cur- 
rent values arc also asymmetric (Fig.[T]3). This asym- 
metry has a different physical origin rooted in its many- 
body excitations, driven by the unequal number of dis- 
crete charge addition and removal channels at opposite 
bias. It is thus clear that the physics of rectification can 
differ widely depending on the strength of the electron- 
electron interaction. 

The non-equilibrium Green's function (NEGF) formal- 
ism is widely established for treating quantum trans- 
port in the SCF regime, for a whole variety of ma- 
terials from nanoscale silicon transistors to molecules, 
nanowircs, nanotubes and spintronic elements. The abil- 
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FIG. 1; (Color online) Experiments showing (a) comparable 
currents reached over unequal voltage widths 0] in the SCF 
limit; (b) unequal currents with comparable widths in the CB 
limit 



ity to incorporate sophisticated quantum chemical mod- 
els d, through averaged potentials makes the NEGF- 
SCF scheme particularly attractive in the community. 
What is not widely appreciated is that this approach does 
not readily translate to the CB regime, even qualitatively 
[13, im, [I4I ■ The CB regime, observed in molecules with 
weak contact coupling [3, Q , manifests clear signatures of 
single-electron charging such as suppressed zero-bias con- 
ductances and abrupt jumps in current. Although several 
approximate treatments [H, [H, [H, [l^ have been sug- 
gested to handle CB within effective one-electron poten- 
tials, the inherent difficulty arises from the fact that the 
open-shell current levels depend on full exclusion statis- 
tics in its many-body Fock space. Even for a minimal 



single-orbital model, it is easy to establish that while the 
open-shell current plateau widths depend on the corre- 
lation strengths, their heights arc independent of corre- 
lation and, in that sense, universal [lO|. The transport 
problem in CB maps onto a rather difficult combinatorial 
problem in many-body space that is hard to capture a- 
priori through a onc-particlc SCF potential, or improve 
upon phcnomcnologically. 

It seems likely that in the limit of weak coupling to 
contacts a proper treatment of the above excitations 
will require solving a set of master equations directly 
in the Fock space of the molecular many-body Hamil- 
tonian [l^, 113] ■ ^ significant penalty is the increased 
computational cost that requires sacrificing the quantum 
chemical sophistication of ab-initio models in lieu of an 
exact treatment of the Coulomb interaction in simpler, 
phenomenological models. Within such an exactly diago- 
nalizablc model, one can capture transport features quite 
novel and unique to the CB regime, such as inelastic co- 
tunneling, gate-modulated current rectification and Pauli 
spin blockade [H, ■ The presence of contact asym- 
metry makes these features even more intriguing, while 
somewhat simplifying the analysis by effectively driving 
the system into equilibrium with the stronger contact. 

In this paper, we first identify the origin of current 
asymmetry with a minimal system of a single spin de- 
generate energy doublet, employing NEGF in the SCF 
limit and master equations for sequential tunneling (26j 
in the CB limit (we refer to this as a "Fock-space" CB 
model). We then extend this "Fock-space model" to a 
general molecular Hamiltonian to explain how multiple 
orbitals in CB allow simultaneous sequential tunneling 
into excited states, making the conductance peak heights 
vary with gate voltage. We also explain the origin of 
exchange in conductance peak asymmetry by the neu- 
tral and singly charged molecule. While this Fock space 
model involves a computationally intensive, exactly di- 
agonalized many-body Hamiltonian about the charge de- 
generacy point arising from the different excitation spec- 
tra accessed, in many cases a simpler approximation 
works. Our analysis of the experimental trends in the 
CB regime reveals that many of the novel signatures such 
as gate dependence of conductance peaks and asymme- 
try flipping can be explained using a simpler "orthodox 
model" that involves an incoherent sum of Fock space ex- 
citations. We conclude by pointing out some limitations 
of such a simplified approach, as well as possible appli- 
cations that may need careful attention to the detailed 
excitation spectrum. 



II. ORIGIN OF CURRENT ASYMMETRIES - 
THE ESSENTIAL PHYSICS 

As mentioned earlier, there are two distinct physical 
limits of transport. In the SCF limit, contact broad- 
enings F arc greater than or comparable with the single 
electron charging U . In the opposite CB limit [/ ^ F and 
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FIG. 2: (Color Online) Origin of Current Asymmetries: a) 
In the SCF regime (stars), when jl S> 7fl for V > 0, the 
left contact keeps the HOMO enoMO level filled to its neu- 
tral charge making conduction possible, while for V <Q, the 
left contact empties the electron leading to the level floating 
out of the bias window, creating unequal plateau widths for 
opposite bias, b) In the CB regime (bold solid line) how- 
ever, for y < 0, charge removal process by the right contact 
is rate determining and occurs in two different ways, while 
for V > 0, charge addition process is rate determining and 
occurs in only one way. This leads to asymmetry in plateau 
heights along the two bias directions, c) While the transi- 
tion between the two limits is hard to accomplish correctly, a 
phenomenological broadening through an artificial enhanced 
temperature (circles) illustrates how the open-shell interme- 
diate CB plateau morphs into a higher effective broadening, 
restoring the SCF result in the limit of F 2> 



single-electron charging dominates. Conductance asym- 
metries in both regimes of transport have been experi- 
mentally observed in molecular conduction. While there 
are ways to handle each regime separately, treatments 
are inherently perturbative, with an approximate treat- 
ment of correlation (in terms of U /F) for the SCF regime, 
and an approximate treatment of broadening (in terms 
of F/C/) for the CB regime. The lack of a small param- 
eter in the intermediate coupling regime {U ^ F) makes 
the exact treatment of transport:, even for a simple model 
system, potentially intractable 

The origin of asymmetric I-Vs can be easily elucidated 
with a minimal model for current conduction through a 
spin degenerate, filled (closed-shell) molecular level dou- 
blet. We assume equal capacitive couplings but unequal 
resistive couplings to the contacts, so that the molec- 
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ular level shifts by half the applied Laplace potential 
and the current onsets arise symmetrically around zero 
bias. In the SCF limit, contact asymmetry results in 
equal currents adiabatically smeared out over a larger 
voltage width along one bias direction than the other. 
This charging based asymmetry has been experimentally 
seen 0, and can be intuitively rationalized as follows. 
Consider a spin degenerate energy level, a highest oc- 
cupied molecular orbital (HOMO) for example, that is 
fully occupied at equilibrium. For asymmetric contact 
couplings 7i ^ 7i?, charge addition dominates for pos- 
itive bias on the right contact and removal for negative 
bias, as shown in Fig. [UJa). For positive bias the en- 
ergy level is maintained at neutrality by the dominant 
left contact and the current flow through the level is de- 
termined by the removal rate. Along reverse bias, in con- 
trast, charge removal by the left contact drives the sys- 
tem away from neutrality towards a net positive charge, 
whose Coulomb cost floats the level out of the bias win- 
dow. This means that a larger bias is needed to fully 
conduct through the level, dragging out the I-V in that 
direction. The direction of the asymmetry flips if conduc- 
tion is through the lowest unoccupied molecular orbital 
(LUMO) instead. Notably the peak currents and initial 
plateau onsets remain the same along both directions, 
but their complete saturations are delayed by different 
amounts. 

The origin and manifestation of current asymmetry is 
qualitatively different in the CB limit, where charge addi- 
tion or removal is abrupt and in integer amounts. Given 
asymmetric contact couplings (7^ 3> 7i?J, the left con- 
tact adds (removes) an electron as soon as the right con- 
tact removes (adds) it, so that the rate determining step 
becomes the dynamics of the weaker right contact. For 
positive bias, charge removal can happen in two ways, 
from ti to t and J,, while for opposite bias the right con- 
tact can add a spin in only one way, either f or | to ti- 
This scheme of charge transfer (Fig. leads to twice 
the current step for positive bias than negative 0, [2^ . 

An important question is whether one can smoothly 
transition from the CB to the SCF asymmetry by pro- 
gressively increasing the broadening. While this is hard 
to do exactly owing to the inherent difficulty involved 
in the broadening many-particle states pl| , for the pur- 
pose of illustration one can add various degrees of ap- 
proximate broadening We choose to do this by in- 
creasing the temperature and incorporating this through 
Boltzmann factors in the many-body occupancies [23| . 
As seen in Fig.[2Ic), this approximate treatment morphs 
the CB asymmetry into the very different version seen for 
the SCF limit. For negative bias on the weaker contact, 
'shell-filling' ^ of the HOMO level with a net positive 
charge creates a CB plateau that is missing in its pos- 
itive bias 'shell-tunneling' counterpart. This extra CB 
plateau gets encased by the broadened manifold, lead- 
ing to the postponed conduction seen in the SCF limit. 
It is worth mentioning though that the correspondence 
is only qualitative and is observed to worsen for higher 



onset voltages, underscoring the inadequacy of thermal 
effects and possibly other phenomenological ways to in- 
corporate broadening in correlated systems, particularly 
in the near equal coupling, non-equilibrium limit which 
combines shell tunneling with shell filling. 



III. COULOMB BLOCKADE FORMALISM: 
FOCK-SPACE VS ORTHODOX 

In this paper, we focus mainly on the CB regime. The 
object of this section is to present the formalism in de- 
tail, beginning with the one that employs the Fock space. 
Here, one needs to keep track not only of various ground 
state charge configurations, but also various excitations 
within each charge state. Following that we derive the 
orthodox model that integrates out the effect of excita- 
tions in an incoherent way. A discussion of merits and 
dc-merits of such a simplification comprises the subse- 
quent sections. 



A. Fock-Space Master Equation 

The starting point is a molecular many-body Hamilto- 
nian 
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where m,p denote the spin-charge basis functions within 
a tight binding formulation, with e, t and U denoting on- 
site, hopping and charging terms, respectively. Exactly 
diagonalizing this Hamiltonian yields a large spectrum 
of closely spaced excitations in every charged molecular 
configuration. The lead molecule exchange processes are 
accounted for within the sequential tunneling approxi- 
mation [23, [2^. This results in a set of master equa- 



tion that involve transition rates i?, 



{N,t)^(N±l,j) 



between 



states differing by a single electron. If one neglects off- 
diagonal coherences, the master equation psj is cast in 
terms of the occupation probabilities Pf^ of each N elec- 
tron many-body state \N,i) with total energy E^^ . The 
end result is a set of independent equations defined by 
the size of the Fock space [2g| 
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For weakly coupled dispcrsionless contacts, parameter- 
ized using bare-electron tunneling rates 7q, [a: left /right 
contact) within a Golden Rule treatment, we define the 
transition-resolved rate constants 



-pNa 



J^\{N~^,J\cJN,^)\' 
7c.|(Af + l,j|ci|iV,z)|2, 



(3) 



4 



where c[^,c^ are the creation/annihilation operators for 
an electron on the left or right molecular end atom cou- 
pled with the corresponding electrode. The transition 
rates are given by 
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riI^/(4'^-M«)- (4) 



for the removal levels {N,i N — and replacing 

(r — > a, / ^ 1 — /) for the addition levels {N, i N + 
fia arc the contact electrochemical potentials, / 
is the corresponding Fermi function, with single particle 
removal and addition transport channels ef^'' = E^^ — 



E^-\ and e 



Na 



El' . Finally, the steady-state 



solution to Eq. ^ is used to get the left terminal current 
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where includes the contributions to R from the left 
contact alone. In the equation above, states correspond- 
ing to a removal of electrons by the left electrode involve 
a negative sign. We usually calculate current in a break- 
junction configuration with equal electrostatic coupling 
with the leads, ^l,_r — Ep ^ eVrf/2. 



B. The Orthodox Model 

In the previous sub-section, the computational com- 
plexity of the master equation defined in Eq. ^ arises 
from the need to keep track of not only charge N , but 
also all configurational degrees of freedom i. Thus the 
evaluation of transition rates defined in Eq. ^ and even- 
tually current (Eq. ^ depends on our knowledge of vari- 
ous many-electron wavcfunctions | A^, i) and total energies 
E^ . The "orthodox" theory of single-electron tunneling, 
however, does not distinguish between excitation levels 
j27l. [28t. Note also that in the orthodox theory, the junc- 
tions are generally denoted by (1,2) instead of {L,R). 
Eq. ([U then becomes 
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We have simplified the second term in the summation 
using a simple change of variables. Let us assume in- 
coherent processes that introduce an additional exclu- 
sion term in the transition process (appendix). Then a 
"golden-rule" calculation gives 
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where D2{E) and Dm{E) are the densities of states of 
the right and middle electrodes, and is the Fermi en- 
ergy of the molecular system. If the level separation be- 
tween one-electron levels is very small, like in a metallic 
quantum dot, the densities of states are approximately 
constant for the calculation of the removal and addi- 
tion rates. By additionally assuming F^" is energy- 
independent, Eq. ([5]) simplifies: 



-AS? 



R,e^ 



1 



exp{AEf/kT)^ 



(9) 



where Rj is the junction resistance. A£^j^ is the tran- 
sition energy corresponding to adding or removing an 
electron at the j*^ contact; from simple electrostatics. 
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where Ci.2,g and Cs are the terminal and total capaci- 
tances respectively, Vu and Vq are drain and gate volt- 
ages respectively, and AU^ is the Coulomb offset for the 
addition or removal of an electron. Notice that we now 
have a simpler analytical and closed form solution to our 
set of master equations, using our trick of summing over 
excitations. In doing so, wc also brought in simpler cir- 
cuit parameters such as junction resistance and capaci- 
tance that aid in our analytical understanding of thresh- 
old voltages and current magnitudes. 

When considering the "orthodox" Coulomb Block- 
ade theory in the regime of strong contact asymmetry 
(i?2 ^ ^i), it is illuminating - and not overly limit- 
ing - to consider what happens at very low tempera- 
tures. Following the analysis of Hanna and Tinkham 
[sot , at low temperatures the ensemble distribution of 
electrons on the middle electrode can be described by a 
delta function where no is the most probable num- 

ber of electrons. The signs of nn and Qq in the follow- 
ing equations will differ from [30| because we are using 
a negative charge carrier convention. The delta func- 
tion probability density reduces Eq. ([7]) to /(Vd,Vg) ~ 
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For low bias the Coulomb 



cost of electrons tunneling across the contacts is high, re- 
sulting in a zero-conductance region limited by the pos- 



itive and negative threshold voltages V^j^ and , re- 
spectively. Outside of this region the transition rates 
simplify: 



I{Vd.Vg) = 
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where sgn denotes the Heaviside sign function. The lin- 
earity of Eq. pip with drain voltage is only interrupted 
when new levels enter into the bias window, causing no to 
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FIG. 3: (Color Online) Illustration of Orthodox Model: The 
orthodox theory parameters define a set of J(V,no) curves, 
represented by the dashed lines. Generally rio is initially zero. 
As a voltage is applied, the current follows the J(V^, no = 0) 
curve until no changes. If |Vcs| < |Vn(,=o|, as is the case 
here, the current will rise linearly out of the zero-conductance 
region; otherwise, there will be a "jump onset." For this 
simulation, Rx = 10 Mf7, 7?2 = 6 Gfi, Ci = 8 aF, C2 = 5.3 
aF, Qo = 0, and T = 2 K. 



change by ±1, which in turn causes the current to "jump" 
in value. 

We thus have an intuitive picture of how I-V (/ — Vd ) 
curves are constructed in the orthodox theory. For given 
system parameters (R, C, etc.) and gate voltage, there 
is a set of /(F, no) curves for different values of no, as 
dictated by Eq. The dashed lines in Fig. [3] corre- 

spond to three members of such a set. In general, both 
no and / are initially zero. As a drain voltage is applied, 
I(V) remains on the no = 0) curve until no changes, 
at which point I(V) jumps to the /(!/, no = ±1) curve. 
Generalizations of Eqs. (5,6a) from Hanna and Tinkham 
[soj allow one to specify the Coulomb Blockade threshold 
voltages, 

y^CB - C^r'(e/2 - r^oe + Qo + CgVg) (12a) 
^CB = Ci\-e/2 - noe + Qo + CgVg).. (12b) 

and the voltages at which the system transitions from 
no to no ± 1 electrons, 

VnP = C^\e/2 + noe - Qo - CgVg) (13a) 
Kb ' = C:2-i(-e/2 + noe - Qo - CgVg). (13b) 

In Eqs. (|12ll3p the positive superscripts refer to positive 
onset voltages, and the negative superscripts are simi- 
larly defined. In Fig. [31 for example, we can see that the 



CB threshold voltage is reached (at 10 mV) before no in- 
creases (at 15 mV), resulting in a linear onset of the cur- 
rent. If, however, were smaller than jVp^'*!, then 
there would be a "jump" onset at the zero-conductance 
region threshold. 

Let us now compare the orthodox and Fock space 
model approaches to transport in the CB regime and ap- 
ply them within the context of experimental trends. 

IV. CB ASYMMETRIES: GATE DEPENDENT 
RECTIFICATION 

One of the simplest consequences of asymmetric con- 
tact coupling is rectification; in other words, a bias direc- 
tion dependence in the I-V characteristics. To calibrate 
with experiments, we not only concern ourselves with rec- 
tification per-se, but also how it is influenced by a gate. 
In fact, experiments [1, 3 showcase gate dependences of 
the rectification properties that are arguably more inter- 
esting than the rectifications themselves. These experi- 
ments (see for example 0, Q) show the following gate- 
able features: (i) a gate-dependent shift of conductance 
peak onsets and (ii) a gate-dependent modulation of the 
corresponding conductance peak heights. In addition, 
there is (iii) a prominent exchange in conductance peak 
asymmetry for gate voltage variations about the charge 
degeneracy point in the stability diagram [3l . We will ar- 
gue that much of the relevant physics has to do with the 
way the molecule accesses various electronic excitations 
under bias, which would require going beyond our one- 
orbital model to a multi-orbital system. Charge addition 
or removal causes jumps in the I-V, while charge redistri- 
bution (excitation) leads to closely spaced plateaus that 
merge onto a linear ramp when summed incoherently. 
In the rest of the section, we will explain how each CB 
model (Fock-space and orthodox) successfully captures 
the gate modulation of the asymmetric I-Vs, as summa- 
rized schematically in Figs. |4] and [T] 

A. Gate Modulation of current onsets and heights 

The onset of conduction is determined by the offset 
between the equilibrium Fermi energy and the first ac- 
cessible transition energy, marked e^'' in Fig. 2] (follow- 
ing the nomenclature in section II A) . This can be varied 
by varying the gate voltage, thereby accounting for the 
variation in conductance gap with gate bias (Fig. ^). 
While the current step and corresponding conductance 
peak are generated by this threshold transition, there fol- 
lows a quasi-ohmic rise in current leading to a subsequent 
constant non-zero conductance in the G-V. This feature 
arises from the sequential access of several closely spaced 
transport channels under bias, arising from excitations 
within the N and — 1 electron subspaces [T^]. While 
net charge addition and removal come at large Coulomb 
prices, excitations involve charge reorganization within 
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FIG. 4: (Color Online) Origin of peak asymmetry and vari- 
ation with gate voltage, (a) For gate voltages that place 
the contact Fermi energy fip m the A'^ electron blockade re- 
gion, the levels align such that fip > ^oo ■ A state transition 
diagram (b) shows the addition and removal of up (down) 
spins resulting in transitions e^"" (bold double arrow) between 
ground states of neutral (light orange) and positively charged 
species (light blue). Also shown in state transition diagram 
are transitions (dashed double arrow) between various config- 
urations of neutral excited state (deep orange) and positively 
charged ground state, labeled e^J" . (c) The resulting G-V 
shows clear asymmetry in conductance peak height due to 
there being more ways to add an electron. Increasing gate 
bias increases the number of excitations available, giving a 
pronounced peak modulation with gate voltage. The inset 
shows the corresponding I-V characteristics [1,11, 

the Fock space that cost much smaller correlation ener- 
gies. 

The presence of multiple orbitals generates several con- 
figurations of excited states, creating more accessible 
transport channels within the bias window. For example, 
in Fig. Ufa) conduction occurs simultaneously via the e^J" 
and e^'' removal channels, ef^'' corresponds to a transi- 
tion between the first excited state '1' of the iV-electron 
neutral species, and the ground state '0' of the N — 1 
electron cationic species. We show four possible configu- 
rations corresponding to the transport channel ej[)'' and 
the corresponding I-V (Fig. HJj) . Increasing the gate bias 
increases both the threshold for current conduction and 
the number of such excited state channels accessed by the 
contacts, thereby altering the height of the corresponding 
conductance peak with gate bias (Fig. HJ;) . 

The previous paragraph illustrates the origin of gate- 



FIG. 5: (Color Online) Experiment from J. Park el al. 0] 
showing gate-rectification properties in the Coulomb Block- 
ade regime: a) Experimental traces, b) Orthodox fit with 
parameters Ci = 0.624 aF, C2 = 0.486 aF, Co = 0.0708 aF, 
Ri = 1 Un, R2 = 75 Un, Qa = -0.05e and T = 2.2 K. c) 
Fit from Fock-space model. 



modulated current as rationalized by the Fock-space CB 
model. One can also explain this within the simpler or- 
thodox models which ignores the identities of the resolved 
excitations by incoherently summing over them. Under 
the approximations of contact asymmetry and low tem- 
perature, the rate is linear in the transition 

energies that increase with drain voltage. With in- 

creasing gate bias one needs a larger corresponding drain 
bias to overcome the zero-conductance regime. At this 
higher drain voltage the coupling has a greater value and, 
consequently, the current magnitude is larger. Physically, 
the drain voltage-dependence of the coupling represents 
a linear approximation of the excitation spectra. Even 
though the orthodox model indiscriminately sums the ex- 
citations within the N and A^ — 1 subspaces, the fact that 
it captures them at all allows it to qualitatively capture 
the modulation of current height. 

Figs.inja) and[6la) show experimental evidence of gate- 
modulation of current onsets and heights. Fig. El^a), from 
J. Park et ai, shows a simple experiment for which the 
negative bias onset is set by moving from the N to N — 1 
electron subspace, while the positive bias onset starts 
where the N electron excitation spectrum moves into the 
bias window. Figs. [nib,c) show the abilities of both the 
orthodox model and the Fock-space model to capture the 
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FIG. 6: (Color Online) Asymmetric CB results showing (a) 
experiment [s^l and (b) orthodox theory with parameters 
Ci = 3.70 aF, C2 = 3.24 aF, Ca = 0.061 aF, i?i = 2 M^, 
R2 = 210 Mfi, Qo = 0.175e, and T = 4.2 K. 



gate-modulated features of the experiment. Fig. [6l^a), by 
Zhitenev et al., is a similar experiment with slightly more 
complex features. Curves (d-h) show the same negative 
bias onset as seen in Fig. [Sj but as the gate voltage is 
further decreased, the negative onset changes to a linear 
onset, representing access to the N electron excitation 
spectrum. The positive bias shows that decreasing gate 
voltage brings the + 1 level closer to the bias window, 
and for curves (a-c) the positive bias type consequently 
becomes a "jump" onset. In spite of these new degrees 
of freedom that must be captured, the orthodox theory 
still models very accurately, as seen in Fig. [DJb). It is 
worth noting that the x-axis in Fig.[B]is Vsd, rather than 
VjjSi which must be accounted for when using Eqs. ^ 
[T3)l of the orthodox theory. A discussion of the extraction 
of orthodox parameters from the experimental curves in 
Figs. O and [5] is included in the appendix. 

Mathematically, it is straightfoward to understand the 
dependence of current height on gate voltage within the 
orthodox theory. Using the experiment of J. Park, et al. 
(Fig. [5^) as an example, we see that there is a jump onset 
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Inserting V^^ ^ for Vd, one electron excited state through transitions et = E^ 



< 



(-) 



CB 

1) curve 



for negative bias voltages. Therefore, 

At V^q \ the I-V transitions to the I{V, nO 
so to find the current height at the onset voltage we can 

use Eq. pT|) for TiQ = ^1- 
finds: 



C1 + C2 
-R2C2C5] 



(Qo + e/2 + CgVg) 



(14) 



Clearly the magnitude of the current at the onset volt- 
age increases with gate voltage (Fig. [SJd), matching the 
experimental result seen in Fig. [5]^ a). 

The accuracy of the orthodox simulation would imply 
at the very least that individual excitations do not play 
an important role in the transport characteristics seen. 
Because the experimental I-V has a strongly linear de- 
pendence on drain voltage seen in Eq. pT|) , it seems that 
the experiments may have measured transport through a 
metallic particle, which has a relatively featureless den- 
sity of states, as is assumed in the orthodox model. 



B. Peak Exchange 

Experiments exhibit a characteristic flipping of con- 
ductance peak asymmetry around the charge degeneracy 
point A in the stability diagram (Figs. [7^, b). Figs.[7]Je,f) 
show typical calculated G-Vs in this regime, featur- 
ing conductance peak asymmetries with respect to volt- 
age bias, arising due to asymmetric contact couplings 
(tl ^ 1r)- Within the Fock-space model, this can 
be explained by enumerating the channels for adding 
and removing electrons under bias (Figs. [7|D,d), with the 
weaker right contact once again setting the rate limiting 
step. The dominant transport channel e^'' corresponds 
to electronic transitions between the neutral and cationic 
ground states [l^, states which SCF theories do take into 
account. In the CB limit, however, there are additional 
electronic excitations that are accessible with very little 
Coulomb cost. These states are responsible for the peak 
asymmetry exchange observed in these experiments, as 
we will now explain. 

The origin of this asymmetry can be understood with a 
simple model system: in our case, a quantum dot with 8 
spin-degenerate levels and = 4 electrons in its ground 
state. When the Fermi energy lies to the immediate right 
of the charge degeneracy point as shown in Fig.[7I^a), only 
transitions between the TV and iV — 1 electron states (4 
and 3) are allowed, with the weaker right contact setting 
the rate-limiting step. For positive bias on the right con- 
tact an electron can be removed from the 4-electron to 
the 3-electron ground state in two ways (Fig. [TJd)). For 
negative bias, however, the electron removed by the left 
contact can be replenished by the right contact back into 
the 4-electron ground state, but also into one of many 
possible excited states e^Q = Ef — E^q , (i > 0) . Since there 
are more ways to bring the electron back (6 shown here), 
the conductance is larger for negative bias (Fig. [TJ;) . The 
situation changes dramatically for a different position of 
the Fermi energy (Fig. [T]:) in the stability diagram lying 
to the left of the charge degeneracy point A with three 
electrons at equilibrium. For negative bias the right con- 
tact adds an electron from the 3 to the 4-electron ground 
state, while for negative bias it returns it to the j*^ 3- 

There now are more ways to remove than add charge 
(Fig. [TJl), so the asymmetry flips (Fig.[7f). 

Analogous to the Fock-space model, the orthodox model 
also captures gate-dependent peak exchange, in spite of 
its approximate treatment of excitations. The origin of 
the asymmetry is once again transitions between the A^ 
to A^ — 1 electron regimes. In Fig.|31 we can see that such 
a change results in a jump onset that has a much higher 
conductance value than a linear onset. Moving from the 
/(V, nO) curve to the /(V, nO — 1) curve, therefore, es- 
sentially captures the excitations of the A^ — 1 electron 
spectrum that are pivotal to the argument in the pre- 
ceding paragraph. The fact that the conductance peak 
switches across zero bias only means that the zero-bias 
state of the system changes from A^ to A^ — 1 electrons. 
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FIG. 7: (Color Online) Origin of exchange in asymmetry of 
conductance peaks: (a), (b) Schematic of vicinity of charge 
degeneracy point A in Coulomb diamond. Also shown are the 
corresponding energy diagrams at threshold. Notice a differ- 
ent set of threshold transport channels between (a) and (b). 
(c) and (d) State transition diagrams illustrating the differ- 
ent excitation spectra accessed on either side of the charge 
degeneracy point A. (e) and (f) G-V plots for scenarios (a) 
and (b) Notice the clear peak-exchange as a result of accessing 
a different excitation spectra in either case. 



which the orthodox method clearly captures. 

Peak exchange has been reported experimentally Q , as 
seen in Fig.[5]^a). Fig.[5]Jb) shows an orthodox simulation 
of the experiment. One can see close qualitative as well 
as quantitative agreement between the experimental data 
and the theoretical simulation. The conductance peaks 
have similar magnitudes, and the exchange of the peak 
asymmetry occurs at —3.75 V in both graphs. The evi- 
dent validity of the orthodox theory in this case demon- 
strates its ability to capture excitation features, as long 
as the features can be linearly approximated. 



V. LIMITATIONS OF THE ORTHODOX MODEL 

In the previous sections, we saw that the orthodox 
model has been fairly successful in reproducing the three 
key trends in asymmetric CB transport; namely, (i) rec- 
tification, (ii) gate-modulation of rectification and (iii) 
exchange of rectification. The agreement between exper- 
iment and theory are quantitatively quite close, indicat- 
ing that the Fock-space model's handling of well-resolved 
excitations is perhaps an overkill, especially considering 
its substantially greater computational complexity. It is 
tempting to conclude that molecules with redox-active 
centers only exhibit incoherent sums of excitations that 
do not manifest well-resolved features. In this section, 
we point out examples where the discrete excitation spec- 
trum can indeed play a noticeable role in molecular trans- 
port experiments, making an orthodox theoretical treat- 
ment quite inadequate. 
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FIG. 8: (Color Online) (a) Experimental trace demonstrat- 
ing peak asymmetry exchange [J]. Compare with Fock-space 
model results in Fig. [3e,f). (b) Orthodox simulation repro- 
ducing peak exchange observed in the experiments. The pa- 
rameters are: Ri = 35 M^, R2 = 350 Mfi, Ci = 0.673 aF, 
C2 = 0.612 aF, Cg = 0.0135 aF, Qo = -0.18e, and T = 4.2 
K. 



Molecules are unique in that they can potentially ex- 
hibit both charge and size quantization Charge 
quantization is enforced by C/ > F, which amounts to a 
large contact resistance compared to the resistance quan- 
tum (the correspondence arises by relating the broaden- 
ing F to an RC time-constant through the uncertainty 
principle, using the capacitance C to set the single- 
electron charging energy U). At the same time, the small 
sizes of molecules make their spectra discrete. While 
transport examples showcasing discrete molecular signa- 
tures are relatively rare, they are pretty commonplace 
in the literature of inorganic quantum dots or 'artificial 
molecules' [s^. Part of the reason is the lower broad- 
ening in these dots, both from the contact as well as in- 
coherent scattering (which is larger in organic molecules 
owing to their conformational flexibility). Recently, a 
novel negative differential resistance (NDR) has been re- 
ported in inorganic double quantum dots [33|, that has 
been explained via the formation of an excited triplet 
state [isl . [s^ l . Other quantum dot experiments routinely 
show Coulomb diamonds with very well resolvable excita- 
tion lines [35I such as due to cotunneling [3^ and Kondo 
correlations [33]. 

Well resolved excitation lines in the conductance spec- 
tra manifest as varying plateaus in the I-V characteris- 
tics. Consider the low bias {Vd < 500 mV) I-V char- 
acteristics of the plot from J-0. Lee et al. [33 , seen in 
Fig. [51 A striking feature is the existence of a plateau at 
onset followed by several other plateaus, sometimes even 
merging into a quasi ohmic rise as a result of several 
unresolvable plateaus. This feature, observed in multi- 
ple experiments from other experimental groups, can be 
easily explained within the Fock-space model by keeping 
track of individual molecular excitations. For example, 
it is well known that the gap between ground and first 
excited states, involving charge addition or removal, is 
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greater than the gap between subsequent excitations in- 
volving charge reorganization. In such a case it is clearly 
seen that a brief plateau occurs at threshold that persists 
until the first excitation is been accessed, as discussed in 
detail in (lo| . 

This leads to an important point. Normally small 
plateaus (few tens of mV) are naturally associated with 
vibronic modes whose energy scales lie in that range. 
Coulomb interactions are usually assumed to be larger in 
energy, because the energy to add or remove an electron 
is much larger. However, the correlation energy to reor- 
ganize charges could actually be very small, are therefore 
quite capable of explaining a range of plateau widths seen 
experimentally, as our calculations show. 

The orthodox theory, on the other hand, cannot even 
qualitatively match the experimental data in Fig. [9l ow- 
ing to its inability to incorporate size quantization effects 
and the associated discrete spectra. From Eq. 1111 it is 
clear that outside of the zero-conductance region and ex- 
cluding jumps due to changes in no, conductance values 
in orthodox theory must remain constant, with a value 
of i?i/i?2Cs. The orthodox theory can capture a plateau, 
and it can also capture a linear rise, but it does not seem 
to capture both in the same I- V curve. Fig. shows 
the best attempt at modeling the experimental data in 
Fig- ma) within orthodox theory; the plateau followed by 
a linear rise seems hard to duplicate. 

A second limitation of the orthodox model comes in 
its treatment of gate voltage. Even though it effectively 
modeled the data in Figs. [5l [H] and [HI there again ex- 
ist experimental features that the orthodox theory can- 
not even qualitatively model. Looking at Eqs. (|12p and 
(fT3|) one can see that a change in gate voltage causes 
a translation in the Coulomb Blockade threshold volt- 
ages; a similar effect, albeit in the opposite direction, is 
seen for the voltage limits at which no changes. Fig. \W\ 
shows how an orthodox I-V curve changes with gate 
voltage for the four possible onset combinations - which 
come from having linear or jump onsets at positive and 
negative bias. The scaling of the conductance gap in 
Figs. Iior a.b) successfully explained the experiments in 
Figs. [5] and [8l However, for experiments with symmetric 
onsets (Figs. [TOh.d). one can see that the orthodox the- 
ory predicts an overall translation in the I-V curve with 
gate voltage. One would think that changing the gate 
voltage would only shift the conducting level closer to 
or further from the bias window, thereby narrowing or 
widening the I-V curve, respectively. Indeed, experimen- 
tal evidence demonstrates such narrowing Q, and the 
Fock-space model captures that quite easily. The ortho- 
dox theory was unable to match this experimental trend 
even qualitatively. 



VI. DISCUSSIONS 

It may seem that the proper treatment of well-resolved 
excitations is rather academic with respect to molecu- 
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FIG. 9: (Color Online) Limitation of Orthodox Model: a) 
Experimental trace [33| showing fine structure, b) The fine 
structure in I-V's, result from keeping track of excitations 
explicitly within the Fock-space CB model, c) An orthodox 
calculation merely maintains the same slope between charge 
addition jumps, thus may not reproduce any fine structure. 



lar experiments. However, we believe that there can 
be more important experimental features that need a 
proper quantitative theory as transport spectroscopy of 
molecules under becomes feasibl e |39ll. In fact, a possible 
explanation for molecular NDRs [40,|41[ could necessitate 
keeping track of excitations in a donor-acceptor molecu- 
lar system [l^; in particular, the notably different life- 
times of ground and excited states, analogous to issues 
relevant to the double quantum dot literature. Small 
molecules could function as tunable quantum dots with 
high single-electron charging energies. Molecular quan- 
tum dots coupled to transistor channels can be impor- 
tant for the detection, characterization and manipula- 
tion of individual spin qubits [39| . the transistor con- 
ductance providing a way to achieve electronic read-out 
[4^ . The large charging energies could allow redox-active 
molecules to operate as storage centers for memory [43| . 
Finally, rectification is important to avoid parasitic path- 
ways in cross-bar logic based architectures [43 |. 

The accurate treatment of well-resolved excitations is 
crucial to molecular operation in this scattering regime. 
But the price paid for this accuracy is the loss of sim- 
plicity associated with orthodox theory. Instead, we 
will need a major improvement in computational algo- 
rithms to handle the exponential scaling of the many- 
body Fock space, such as partial configuration interac- 
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FIG. 10: (Color Online) Gate-voltage dependence of I-V 
curves in the orthodox theory. Blue arrows indicate the direc- 
tion of movement of the onset voltages with increasing gate 
voltage. The four plots correspond to the four possible onset 
combinations: (a) a jump onset at negative bias and a lin- 
ear onset at positive bias, for which the I-V conductance gap 
widens with increasing gate voltage; (b) a linear onset then 
a jump onset, with the gap narrowing around zero-bias; (c) 
two linear onsets, which corresponds to a translation of the 
I-V curve, and (d) two jump onsets, which correspond to a 
translation in the opposite direction of (c). 



tion (CI) with frontier orbitals coupled to the standard 
master equations. Major challenges involve the proper 
inclusion of the contact boundary conditions, and doing 
justice to broadening [2l|. Also, one could incoherently 
sum selected excitations and resolve the relevant ones, 
if we have some prior idea which ones could be impor- 
tant. Spin unrestricted approaches may capture some of 
the salient effects, although they are unlikely to capture 
the current heights ( which depend on difficult parameter- 
independent combinatorial arguments), and even more 
significantly, the number of current plateaus (since an 
energy-independent one-electron potential does not gen- 
erate enough poles). Needless to say, there is enormous 
room for theoretical activity and device potential in this 
domain. 

It is worth emphasizing that the wide success and 
popularity of the SCF-NEGF approach often leads re- 
searchers to automatically assume that the method will 
work in the CB regime [H, El]. The NEGF method 
is indeed prolific in its flexibility to incorporate contact 
microstructurc, molecular chemistry, broadening, interef- 



erence and electrostatic details. However, the method, 
as widely implemented, is fundamentally limited by the 
need to work with a one-electron potential, at best in- 
corporating many-body effects as corrections to this po- 
tential. While this may be sufficient to capture equi- 
librium properties like total energy, transport measure- 
ments can potentially probe the rich excitation spectra 
that are hard to incorporate readily into the average elec- 
tronic potential, even parametrically. For a typical self- 
energy, there simply are not enough poles in the one- 
electron Green's function to do justice to the many more 
many-electron transitions characterizing transport, not 
to mention the complications of the full counting statis- 
tics in Fock space under non-equilibrium that show up 
as effective interaction-independent degeneracies of the 
open-shell levels. 

In this paper, we have discussed contact asymmetry 
as a paradigm to delve into various experimental fea- 
tures, outlining the qualitative difference and cross-over 
between the weakly correlated SCF and the strongly cor- 
related CB regimes. An understanding of the asymme- 
try was provided using a basic model, and then extended 
to conjugated molecular systems. Novel experimental- 
trends were identified in the CB-regimc of transport. 
Two different approaches, the Fock-space CB and the 
orthodox CB, were discussed in detail. While a simpler 
orthodox approach captures the salient asymmetric ef- 
fects, the Fock-space approach may be essential in ex- 
amples where the interplay between specific excitations 
governs the principle operating physics and possibly also 
interesting device applications envisaged. 
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VIII. APPENDIX I: EXTRACTING 
PARAMETERS FOR ORTHODOX MODEL 

The low temperature, strong contact asymmetry ap- 
proximation to the orthodox theory, as introduced in sec- 
tion III; is most useful because it simplifies the process 
of simulating experimental data. Whereas the general 
orthodox formalism requires best-fit computational tech- 
niques, Eqs. pn -12) allow for a straightforward analytic 
calculation of the system circuit parameters. Consider 
Fig. (O as an example of data to be simulated with the 
orthodox theory. To start, consider the Vg = — 0.85F 
I-V curve (light blue). We have a linear onset at ap- 
proximately 19 mV , and a "jump" onset at, say, —24 
mV . Plugging these values into Eqs. (11a) and (12b), 
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respectively, and solving finds: 



Ci{l9mV) = C2(24mF), and 
' + 2Qo + 2CgVg 



2(19771^) 



(15) 
(16) 



Consideration of a second I-V curve in the data gives us 
further information; here we will use the Vq = —0.55^ 
(green) curve, which has a positive onset voltage of ap- 
proximately 53 mV. After plugging the correct onset and 
gate voltage values into Eq. (11a) for both I-V curves 
considered, one can solve for a new, independent equa- 
tion: 



AVg 



(53 - 19)mV 



Cq — 8.82Cq 



(17) 



We now have three equations and four unknowns (Ci, 
^^2, Cg, and Qq). There are now two approaches one 
can take. The first, which was applied in the simulation 
shown in Fig. 4(b), is to simply take a small, reasonable 



value of Qoj such as — 0.05e, and solve the remaining 
equations for the other variables. However, this will not 
always work, for a reason that is somewhat subtle. It is 
important to note in Fig. 4(a) that not only does the 
Vg = — 0.85y I-V curve have a linear onset at 19 mV , 
but it also does not have a jump onset at positive bias 
until at least 100 mV . Setting Eq. (12a) greater than 
100 mV provides a fourth restriction. Indeed, for the 
experiment in 4(a) this is trivially satisfied; Fig. [6l on 
the other hand, is an example of data for which such a 
fourth equation is necessary for a successful simulation. 

A final parameter one can obtain is i?2. From Eq. ([7]) 
we know that the conductances of the linear portions of 
orthodox curves have the value CxjRiC^. Having al- 
ready calculated all of the capacitance values it is thus 
straightforward to find the slope of the experimental I-V 
curve and find i?2- cannot be found so easily, but its 
relative unimportance in the regime of contact asymme- 
try means that a simple estimate of R\ ~ i?2/100 will 
usually suffice. 
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